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Abstract 

This paper deals with different ways to extract the effective two-dimensional lower level dynamics 
of a lambda system excited by off-resonant laser beams. We present a commonly used procedure for 
elimination of the upper level, and we show that it may lead to ambiguous results. To overcome 
this problem and better understand the applicability conditions of this scheme, we review two 
rigorous methods which allow us both to derive an unambiguous effective two-level Hamiltonian of 
the system and to quantify the accuracy of the approximation achieved: the first one relies on the 
exact solution of the Schrodinger equation, while the second one resorts to the Green's function 
formalism and the Feshbach projection operator technique. 
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I. INTRODUCTION 



When dealing with complicated multilevel systems, such as atoms or molecules, it is 
necessary to look for allowed restrictions of the Hilbert space which can lead to simplifications 
of the computational work. Thus, one usually forgets states which are not populated initially 
and not coupled, either directly or indirectly, to initially occupied states. 

It is sometimes possible to go beyond this first step and isolate some subset of initially 
occupied states if they are only weakly and non-resonantly coupled to the others. The 
effective dynamics of such a subset can then be approximately described by a Hamiltonian 
of smaller dimensions than the original one, in which the effect of couplings outside the 
relevant subspace is accounted for by additive energy shifts and couplings. The procedure 
which allows one to get rid of the irrelevant states and derive this effective Hamiltonian is 
called adiabatic elimination. One of the simplest examples of such a situation is the case 
of a three-level lambda system, the low levels of which are initially populated and non- 
resonantly coupled to the initially empty upper level via detuned harmonic perturbations: 
through adiabatic elimination, it is possible to reduce the problem to an oscillating two-level 
system, which has been widely studied, for instance in atomic physics 

The aim of the present paper is to understand how the adiabatic elimination procedure 
works on this simple example and how it should be performed. After briefly presenting the 
model (Sec. |H}, we show that an Ansatz commonly used in the literature to adiabatically 
eliminate the excited state may lead to ambiguous results (Sec. IIIIJ) . In order to better 
understand this scheme and explicit its conditions of applicability, we then review two rig- 
orous approximation methods to treat the problem: the first one relies on the solution of 
the Schrodinger equation and consists in neglecting the fast oscillating terms in the exact 
expression of the amplitude of the excited state, which is then injected back into the dynam- 
ical equations for the lower states (Sec. II V|) : the second one resorts to the Green's function 
formalism and makes use of the pole approximation which is discussed in detail (Sec. IV|) . 



II. THE LAMBDA SYSTEM 

In this paper, we shall consider an atomic lambda system[9j] consisting of two lower 
states {\a) , \b)} coupled to an excited level |e) via two off-resonance lasers: the detunings 
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FIG. 1: Lambda system's level scheme: a) in the Schrodinger picture, b) in the Rotating frame. 



are denoted 5 k = uj k — (k = a, b) where frw k = E e — E k and ui^ are the frequencies of 
the lasers (cf figure QJi). 

Turning to the rotating frame defined by the transformation — ► \(p) — e +i & \ , 
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where £ = 







10 a - L0 b - 
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where A = Sa + Sb ; and (f2 a ,f2 b ) denote the Rabi frequencies of the lasers coupling |a) to \e) 
and 1 6) to |e), respectively. Note that we have implicitly chosen the origin of the energies 
between the two states \a) and \b) (see figure GJd). 

From now on, we shall assume |A| 3> \8\ , \Q a \ , |0&|. If the system is initially prepared in 
a superposition a \a) + j3 \b), the excited state will then essentially remain unpopulated, 
while second-order transitions will take place between the two lower states: this constitutes 
the so-called Raman transitions, which play an important role in atomic and molecular 
spectroscopy, and have recently become importantprocesses in laser cooling and trapping 
[2] and in quantum computing proposals with ions 3], atoms [4j, and solid-state systems jj 
In these conditions, it is natural to restrict the Hilbert space to the relevant states \a) and 
I b) and to describe their dynamics by a 2 x 2 effective Hamiltonian H e ff. In the following 
sections, we describe and discuss different ways to derive H e ff. 



III. ROUGH ADIABATIC ELIMINATION IN A LAMBDA SYSTEM 



The usual way to eliminate the excited state from the Schrodinger equation, written for 

a 



the state 
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ia (t) 

z/3(t) = f/8 + ^7 

z 7 (t) = + f(3 + A 7 



(2) 



consists in claiming 7 (t) = 0, which implies, by solving the last equation 

1 = tCh rP 3 

and injecting (j^J back into the dynamical equations for a and (3, which yields 



ihdt 



where the effective two-level Hamiltonian writes 
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If we had performed the same calculation in a shifted picture defined by the transforma- 

we would have obtained 
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and, applying the same Ansatz 7^ (£) = as before (we assume 77 7^ —1), we would have 
derived 

Mb 
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which yields the effective Hamiltonian H e ff ;7) in the shifted picture. Finally, subtracting 77 A 
from H e ff v we would have obtained the effective Hamiltonian H e ff jT) 
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We are thus led to the obviously unphysical conclusion that the effective Hamiltonian, 
i.e. the effective dynamics of the system depends on the picture where the elimination is 
performed. This raises questions about the Ansatz we used: in which picture, if any, does 
it apply, and what is its physical meaning ? To understand better when and how to employ 
this scheme, we investigate two rigorous elimination methods in the next sections, which 
yield both an unambiguous expression of the effective Hamiltonian and the level of accuracy 
of the approximation achieved. 



IV. RIGOROUS ADIABATIC ELIMINATION IN A LAMBDA SYSTEM 
THROUGH SOLUTION OF THE SCHRODINGER EQUATION FOR THE AM- 
PLITUDE OF THE EXCITED STATE 

In this section, we propose a straightforward elimination scheme based on the analysis of 
the exact expression of the amplitude 7 of the excited state: resorting to a simple mathemat- 
ical argument, we identify its relevant part 7 re z which mainly contributes to the dynamics 
of (a, f3); injecting it back into the dynamical equations, we then derive an unambiguous 
expression for H e ff which yields an approximation of the dynamics of the system, the accu- 
racy of which can be quantified. This method moreover allows us to specify the applicability 
conditions of the previous scheme. 



A. Exact solution of the problem 

The derivation of the exact solutions of the Schrodinger equation can be straightforwardly 
performed by finding the eigenenergies of the system and using the boundary condition 
\ip (t — 0)) = a I a) + Pq \b) to determine the coefficients of the Fourier decompositions of 
the different amplitudes. We do not reproduce these calculations but only summarize the 
results which are useful for our purpose. Introducing the reduced variables (A, Xk= a ,b) such 
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TABLE I: Expansions of the different parameters of the exact solutions. 
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that Ae = \ k=ab e = £|, with < e < 1, (A, Afc =a ,fe) = O (1), and A > 0[1C|, one readily 
shows 

' «(*) = ELi^ e ~ iAa:fet 

< P(t) = ELiB k e-^ (8) 
7 (t) = ELi C k e~^ 

v 

where {x/c} fc=1 2 3 are the solutions of the equation 

x 3 - x 2 - (A 2 + |A a | 2 + |A b | 2 ) e 2 x + A 2 e 2 + Ae 3 (|A a | 2 - |A 6 | 2 ) = (9) 

and the coefficients {A k , B k , C k } k=1 2 3 are determined by the boundary conditions. Table 
IJ displays the expansions in e of these different parameters. The rigorous solutions are 
represented in figure in the regime |A| ^> \S\ , |O a | , the amplitudes a and (3 show 
an oscillating behavior, which, as expected, is much alike a two-level Rabi oscillation; 7 
oscillates much faster than a and (3. The comparison between the exact solutions and the 
approximations obtained in the previous subsection (see figure OJ) shows that the previous 
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FIG. 2: Exact solutions for a, [5 and 7. The values of the parameters used for the simulation are: 
a = V /I73, fa = ^273, 6/ A = 0.1, Qi/A = 0.1 x e~^/ 3 , 2 /A = 0.1 x e^/ 2 . 

elimination method provides a valuable approximation only in what we called "the natural 
picture" . In the following, we clarify why this is so through a straightforward and rigorous 
elimination scheme. 



B. Straightforward adiabatic elimination 

1. Preliminary remark 

Consider the differential equation 



if{t)± 5 -f{t) = ^{Ke-^ + A^) 



(10) 



the solution of which writes 
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FIG. 3: Comparison of the exact solutions with the rough approximation calculated in a shifted 
picture, for different values of the shift 77. The values of the parameters used for the simulation 
are: a = y/T/3, (3 = y/2/3, 6/ A = 0.1, Ui/A = 0.1 x e^/ 3 , U 2 /A = 0.1 x e~ i7r / 2 ; exact solution 
(full line), 77 = (dashed line), 77 = 0.3 (broken line), 77 = 3 (longbroken line). 

If one is interested in a solution valid to first order in e, the O (e 2 ) term can be neglected, 
or, equivalently, one can forget the corresponding term directly in the differential equation 
which then reduces to 

l f(t)± 6 -f(t) = ^e-^. 
In the same way, one straightforwardly shows that the solutions of the equation 

if (t) + (±i - rjA^j f (t) = I (A.e-^X + A A e-^ 1+ ^) (11) 

where 77 is an arbitrary real number and all the other parameters obey the same relations 
as before, coincide (to order O (e)) with those of the simpler equation 

if (t) + (± 5 - - rjA^j f (t) = ^e-^^. 
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2. Rigorous adiabatic elimination of the excited state 

Let us now return to our initial problem. Injecting the exact expression (jHJ) of 7 into (j2J), 
we get the equations 

ia = Aa + cJ^e-* XlA * + C 2 %~ iX * A 1 + C 3 %-^ At 
2 2 2 2 

which are of the form (jlUJ) . with 

Neglecting the last term of each of these equations, or equivalently, replacing 7 by its "rele- 
vant" component 

2 

lrd (t) = 7 (*) - C 3 e-^ At = ^ C fc e-^ A * 

fc=i 

then leads to approximate forms for (a, valid to order O (e) terms. Let us now relate 7 re / 
to a, /3: we have 
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we deduce from (|T2j) that 

Jrd = -^*-^ + 0(<?). (13) 

This expression coincides with (J3J): it means that, in the natural picture, the rough adiabatic 
elimination procedure indeed leads to a valid approximation of the actual dynamics of the 
system, to order O (e), which agrees with the remark we made about figure 01 and legitimates 
the expression (JHJ) for the effective Hamiltonian. 

If we turn to the shifted picture defined by the transformation \cp) — > \lp v ) = e~ ir?A * \cp), 
the exact expression of the amplitude of the excited state is now 

3 
k=l 

and leads to the dynamical equations 

V 2 / 2 2 2 

which are of the form (fTTjl . By similar calculations as above, one shows that a good approx- 
imation of the dynamics of the system is then obtained through replacing 7 by its relevant 
part 

^ el = -^a-^=e-^ re , (14) 

This expression differs from the rough adiabatic elimination result ((Zj). As can be easily 
checked, however, when rj = O (e), ((7|) and (|14|) only differ by O (e 2 ) terms which are not 
significant at the level of accuracy we consider. The rough elimination procedure thus 
appears as a practical (though physically not motivated) trick which works as the long as 
the origin of the energies lies in the " neighbourhood" of the midpoint energy i?o of the two 
lower states, or, to be more explicit, as long as E / A = O (e). 

The rigourous method we have employed here allowed us to clarify the applicability 
conditions of the rough procedure but it requires the exact solution of the Schrodinger 
equation. In the next section, we present a more systematic and elegant approach based on 
the Green's function formalism which can be used to generalise the results presented here 



to more complicated level schemes 



8]. 
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V. RIGOROUS ADIABATIC ELIMINATION IN THE GREEN'S FUNCTION 
FORMALISM 



A. Overview of the Green's function formalism 



This overview summarizes the basic features of the Green's function and projection op- 

fl 

erator formalism first introduced in [y] and extensively presented in |7[. Given a system of 
Hamiltonian H = H + V, comprising a leading part H and a perturbation V, one defines 
the Green's function G (z) as follows 
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The evolution operator of the system can then be derived through the formula 
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where C + and C_ denote two parallel lines just above and below the real axis, oriented from 
the right to the left and from the left to the right, respectively. 

Let V be the subspace spanned by some relevant eigenstates of H , and let P and Q = 
I — P be the orthogonal projectors on V and V ± = Q, respectively. Then, one can show 

PG (z) P P 



z - PH P - PR (z) P 
where the displacement operator R (z) is defined by 

Q 



R(z) = V + V —QfJ a Q _ QVQ 

B. Application to the lambda system 

In our case, 

P = \a) (a\ + \b) (b\, Q = |e) (e| , H = h 



V. 



s 

2 
















"a 

















5 
2 





2 
















A 






n b 






whence 



PG (z) P 



P 



M z) 



M(z) = z + XehAa z - 



(HAe) 2 
z-HA 



|A a | A*A 
A a A^ | Aft 



2 



11 



To compute 



PU(t)P=— [ dze- izt/h P/M 
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(15) 



one has to calculate the residues of e tzt / h P/M (z) at its poles which are the solutions of 
det [M (z)) = 0. Setting x = one gets 
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which coincides with © and thus leads to the same results {xk} k=1 2 3 summarized in Table 
HI The associated residues {TZk} k=1 23 are readily found to be 
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In both cases, the third pole only contributes to the second order in e to the evolution 
operator (fTHjl . whereas the first two are O (1). To the first order, the last pole can thus be 
omitted: this constitutes the so-called pole approximation. 

A simple and straightforward way to implement the pole approximation is to replace 

M (z) by 



z + XehAcr z + 



jhAe) 

hA 



|A a | A*Afc 
A a A^ \Xt\ 



in (|15|) which boils down to replacing PR (z) P by PR (0) P in the expression of PG (z) P. 
One readily shows that, as desired, the replacement of PR (z) P by PR (0) P in PG (z) P 
discards the irrelevant pole and its associated residue, while leaving the others unchanged 
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(up to (e 2 ) terms). Moreover PG (z) P can now be put under the form PG (z) P ~ ^ 
where H e ff = PH P + PR (0) P is a 2 x 2 Hermitian matrix, independent of 2 which can 
be interpreted as the effective Hamiltonian governing the dynamics of the reduced two- 
dimensional system {\a) , \b)}: a simple calculation again leads to the expression (JHJ) for 

H eff- 

Let us finish this section by some remarks. The calculations above have been performed 
assuming that the energy E at the midpoint of the two ground states spanning V is zero: 
if one shifts the energies so that E ^ 0, all the previous expressions and calculations hold, 
up to the replacement of z by its translated z + E , as can be easily checked. The expression 
of the effective Hamiltonian thus becomes 

H eff = PH P + PR(E )P, 

which is of course consistent with the particular case (E = 0) considered above. It is in- 
teresting to note that H e ff always takes the same form, wherever one chooses the origin of 
the energies, i.e. the sum of the projected unperturbed Hamiltonian and the operator PRP 
evaluated at the middle energy of the subspace V. Finally, let us also note that if PRP 
is evaluated at an energy E\ = E + O (e), the previous approximation remains valid, as 
the residues will only be affected by O (e 2 ) terms, which are not significant at the level of 
accuracy considered. 



VI. CONCLUSION 



The goal of this article was to clarify the scheme usually employed to derive the effective 
two-dimensional Hamiltonian of a lambda system excited by off-resonant lasers: in particular 
we have reviewed two methods which enabled us to rigorously derive the effective dynamics 
of the system, up to terms to well-established order of magnitude in a small parameter; 
this study also allowed us to specify the applicability conditions of the rough elimination 
procedure. The second of these schemes relies on the Green's function formalism and the use 
of projectors on the different relevant subspaces of the state space of the system: this fairly 
general and elegant tool naturally leads to generalisations to more complicated multilevel 
systems, as shall be considered in a forthcoming paper. 
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